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GTPases are molecular switches that regulate a wide-range of cellular processes. The GPN-loop GTPase (GPN) is a sub-family 
of P-loop NTPase that evolved from a single gene copy in archaea to triplicate paralog genes in eukaryotes, each having a 
non-redundant essential function in cell. In Saccharomyces cerevisiae, yGPNI and yGPN2 are involved in sister chromatid 
cohesion mechanism, whereas nothing is known regarding yGPN3 function. Previous high-throughput experiments 
suggested that GPN paralogs interaction may occur. In this work, GPN|GPN contact was analyzed in details using TAP- 
Tag approach, yeast two-hybrid assay, in silico energy computation and site-directed mutagenesis of a conserved Glu 
residue located at the center of the interaction interface. It is demonstrated that this residue is essential for cell viability. 
A chromatid cohesion assay revealed that, like yGPNI and yGPN2, yGPN3 also plays a role in sister chromatid cohesion. 
These results suggest that all three GPN proteins act at the molecular level in sister chromatid cohesion mechanism as 
a GPN|GPN complex reminiscent of the homodimeric structure of PAB0955, an archaea! member of GPN-loop GTPase. 



Introduction 

Cohesins play a central role in pairing of duplicated chromo- 
somes''^ and in a wide-range of other functions, such as double- 
strand DNA repair** and transcriptional regulation.*'^ For efficient 
activity of cohesins, numerous additional proteins are required. 
Recently, two proteins, namely the GPN-loop GTPases GPNl 
and GPN2, have been shown to be involved in cohesion dissocia- 
tion^ and establishment,' respectively. 

GPN-loop GTPases are found in all sequenced eukaryotes 
and in almost all archaea but not in bacteria.'" In eukaryotes, 
three distinct homologs (GPNl, GPN2 and GPN3) are found, 
whereas a single copy is found in archaea. This family of con- 
served GTPases has been named GPN-loop GTPases (GPN) 
due to the presence of a conserved Gly-Pro-Asn motif inserted 
into the GTPase core-fold." These GPN-loop GTPases belongs 
to monophyletic GTPases superclass within P-loop containing 
nucleoside triphosphate hydrolases fold, which is divided into 
two large groups designated TRAFAC (after translation factors) 
and SIMIBI (after signal recognition particle, MinD and BioD) 
(see review in ref 12). This later group includes, among others, a 
group of metabolic enzymes with kinase activity, signal recogni- 
tion particle (SRP) GTPases and the assemblage of MinD-like 
ATPases, which are involved in protein localization, chromosome 
partitioning and membrane transport.'^ Up to now, the three- 
dimensional structure of a GPN-loop GTPase, that of PAB0955 
from the hyperthermophilic archaeon Pyrococcus abyssi, has been 
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determined. Its structural study revealed a general fold related 
to SIMIBI GTPase." The crystal structure also shows that this 
archaeal GPN is homodimeric." In the yeast, Saccharomyces cere- 
visiae, individual deletions of either gene, YJR072C (yGPNI), 
YOR262W (yGPN2) or YLR243W (yGPN3), are lethal.'" Thus, 
these three paralogous GTPases fulfill an essential function and 
are not functionally interchangeable. 

The present work sheds new lights on yGPNs function both 
at the molecular and cellular level. Pull-down experiments using 
yGPNI as bait first showed that in yeast yGPNI can be associ- 
ated with yGPN2 or yGPN3. Then, in order to characterize the 
molecular basis of these complexes, and in the absence of crys- 
tallographic data, comparative models of the three yGPN para- 
logs were built using PAB0955 as a template, and modeling of 
yGPNl|yGPN2 and yGPNl|yGPN3 complexes were performed. 
These complexes revealed putative critical residue at the interface 
of two monomers: Glull2 for yGPNI and yGPN2, and GlullO 
for yGPN3. Mutations of these residues were realized, and it 
was demonstrated using a two-hybrid system assay that they are 
essential for paralogous interactions. Based on a sister chromatid 
cohesion assay, we also found that overexpression of yGPN3 pro- 
tein promotes sister chromatid separation during anaphase. 

Results 

yGPN2 and yGPN3 belongs to the yGPNI interacting net- 
work. Large-scale two-hybrid experiments in S. cerevisiae 
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Figure 1. Comparative model of theyGPNI (red)|yGPN3 (blue) interac- 
tion complex. Secondary structure elements are indicated by arrows for 
p-strands and helicoidal springs for a-helices. Top panels show close-up 
views on the conserved Glu residues located at the yGPNI |yGPN2 (left) 
and yGPNI |yGPN3 (right) interfaces. The Ca trace of molecules is shown 
as tubes, whereas amino acid side-chains are shown using colored 
balls and sticks (black for carbon, red for oxygen and blue for nitrogen). 
Hydrogen bonds made between the conserved Glu residues of one 
partner with backbone atoms of the other partner are shown as gray 
dots. Each Glu residue makes three side-chain-to-main-chain hydrogen 
bonds except in the yGPNI |yGPN2 dimer where a proline is present in 
yGPN2 (no nitrogen backbone atom available). The transparent mo- 
lecular surface was built using the surf tool of VMD,^' whereas the top 
panels were built using Molscript" and rendered using RasterSD.^^ 



suggested that yGPN2 and yGPN3 proteins could interact with 
yGPNI."'"' Recently, the co-purification of tagged yGPNI with 
yGPN2 has been demonstrated by mass spectrometry in S. cerevi- 
siaeF Here, we show that yGPNI |yGPN3 interaction could also 
be detected using tandem affinity tag strategy coupled to a shot- 
gun mass spectrometry approach. For this, yGPNI TAP-tagged 
protein was chosen, and fiarmaldehyde was applied on intact cells 
harvested at the exponential phase to freeze transient interactions. 
Purified complexes were subjected to a short-run electrophoresis 
in denaturing condition and then in-gel proteolyzed with tryp- 
sin. The resulting peptides were identified by nanoLC-MS/MS. 
A total of 815 different peptides were detected. Table 1 reports all 
the proteins identified with at least two different detected pep- 
tides and filtered of contaminants frequently observed with the 
TAP-TAG assay. Among these, 16 proteins are GPN proteins as 



well as proteins belonging to the transcription machinery, six of 
which are also found in conditions without the cross-linker folm- 
aldehyde (Table 1). As expected, the yGPNI is the most detected 
protein, with about 15 different peptides observed by nano 
LC-MS/MS. Besides RNAPII subunits, Rpbl, Rpb2, Rpb3 and 
Rpb5, and yGPN2 protein, which have already been observed 
as interacting with yGPNI in previous TAP-Tag studies,'^ 10 
new interacting proteins were found co-purified with yGPNI: 
two subunits of RNAPI (Rpal90 and Rpal35), one subunit of 
RNAPII (Rpbll), two subunits of RNAPIII (Retl and Rpo31), 
two RNA polymerase subunits common to RNA polymerases I 
and III (Rpc40) and RNA polymerases I, II and III (Rpb8), RTRl 
(Rpbl phosphatase), Rba50 (RNA polymerase Il-associated pro- 
tein) and yGPN3. We observed that most of the GPNs interac- 
tions were only detected in the presence of a cross-linker, which 
may allow the freezing of transient or weak interactions. 

Comparative modeling of yGPNI, yGPN2 and yGPN3. 
To understand the molecular basis of possible direct interac- 
tions between yGPNI |yGPN2 and yGPNI |yGPN3, molecular 
models were sought for each of these proteins as well as their 
putative complexes using our in-house modeling pipeline. The 
closest identified protein template in the Protein Data Bank"* is 
the structure of the archaeal protein PAB0955 from Pyrococcus 
abyssiP The sequence identity between yGPNs and PAB0955 is 
relatively low, around 20%, but there is a significant similarity in 
the first half of the sequence alignment, which gave us high confi- 
dence in the choice of the template. Stereochemical and energetic 
evaluations of the three molecular models are shown in Table 2. 
Values indicate that all three models are of good quality compared 
with the 2.3A-resolution crystallographic structure of PAB0955. 
The template PAB0955 was described as being a homodimer," 
and we wondered whether heterodimers yGPNl|yGPN2 and 
yGPNl|yGPN3 are plausible. For this, we built molecular mod- 
els for each of these heterodimers as described in the methods 
section. This stage was straightforwardly executed, taking advan- 
tage of the dimeric structure of the template PAB0955. Figure 1 
displays the yGPNI |yGPN3 complex. Very few adjustments were 
required to fit side chain conformations in both heterodimers. 
Root mean-square deviations of Ca atoms (Ca-RMSDs) were 
calculated after a global superimposition of these complexes onto 
the dimeric PAB0955 using the sup3d software program,^" and 
both gave 0.26 A for yGPNl|yGPN2 and yGPNl|yGPN3 dimers. 
The computed buried surface area of yGPNl|yGPN2 is 3860 A" 
and 3706 A^ for yGPNI |yGPN3 compared with 3988 A^ in the 
dimeric PAB0955. All these parameters indicate that heterodi- 
mers yGPNl|yGPN2 and yGPNl|yGPN3 present energetic and 
geometrical values close to that of the dimeric template. 

The conserved residues E112 in yGPN2 and EllO in yGPN3 
are critical for in vivo yGPNl|yGPN2 and yGPNl|yGPN3 
interactions. A detailed analysis of the dimeric structure of the 
archaeal PAB0955 protein revealed the presence of a buried 
charged residue at the dimeric interface, the glutamate at posi- 
tion 107 (El 07), which corresponds to El 12 in yGPNI, El 12 in 
yGPN2 and EllO in yGPN3. This amino acid is located right 
next to the G3 box belonging to the core structure of the GTPase 
active site. It is strictly conserved among the GPN-GTPase 



464 



Cell Cycle 



Volume 12 Issue 3 



Table 1. Spectral count differences detected with or without cross-linker 



Protein 


Accession 


name 


Number 


yGPN1 


YJR072C 


yGPN2 


YOR262W 


yGPN3 


YLR243W 


RPA190 


YOR341W 


RPA135 


YPR010C 


RP021 


YDL140C 


RPB2 


Y0R151C 


RPB3 


YIL021W 


RPB11 


YOL005C 


RET1 


YOR207C 


RP031 


Y0R116C 


RPB5 


YBR154C 


RPC40 


YPR110C 


RPB8 


YOR224C 


RBA50 


YDR527W 


RTR1 


YER139C 



witiiout 
crcsslinlter 
exp 1 


% 

coverage 


No. of {*) 
peptides 


Score 


33 


13 


659 


- 

8 
12 

9 


- 

12 
11 
2 


- 

502 
414 

78 


5 


1 


34 


4 


1 


23 





without 






crosslinl<er 






exp 2 




% 

coverage 


No. of 
peptides 


Score 


33 
- 


12 

- 


605 
- 


4 


7 


216 


4 


3 


94 


9 


2 


86 


5 


1 


34 









witii 
crosslinlter 
0.15% 


% 

coverage 


Nn nf 

peptides 


Score 


46 


18 


971 


11 


3 


156 


21 


3 


169 


7 


10 


403 


9 


8 


342 


18 


26 


1322 


34 


37 


1600 


14 


3 


190 


2 


2 


55 


13 


2 


151 


3 


1 


38 


7 


3 


132 











with 






crosslinker 






0.20% 




% 

coverage 


Nn nf 

peptides 


Score 


37 


19 


1199 


34 


8 


558 


27 


5 


343 


28 


34 


1898 


24 


22 


1114 


38 


51 


3077 


50 


58 


3341 


44 


10 


627 


50 


4 


186 


15 


12 


543 


2 


2 


92 


33 


5 


353 


25 


6 


335 


17 


2 


119 


13 


6 


371 


33 


5 


352 



I , subunit of RNA POLI; , subunit of RNA POLII; ■, subunit of RNA POLIII; ■, subunit commun to RNA POLI, II and III; ■, other 



(Fig. SI). The presence of several hydrogen bonds between 
side-chain atoms of GIu residue and several nitrogen backbone 
atoms suggests that this residue plays an important role in the 
energetics of the complex interface (Fig. 1). Preliminary in silico 
mutagenesis was performed at this position with every other 19 
amino acids (AA) using in-house software.^' It was found that 
most AA replacements were deleterious. The energy difference 
AE (E^^^ - E^^) for lysine (K) replacements in yGPNl|yGPN3 
complexes are indicated in Table 3. It shows that the double 
mutants yGPNl,^„,^)|yGPN2,^„,^), yGPNl,^jj,^JyGPN3(£„„K) 
are clearly less stable than single mutants. It should be noted that 
in their respective monomeric form, yGPNl^,|^jj2' yGPN2^,|^jj2 
and yGPN3j,,^|ug are fully solvent-exposed. Consequently, com- 
putational replacement of these residues in monomeric GPN pro- 
teins did not reveal disruptive effects. To experimentally confirm 
that GlullO/112 is critical for GPN complexes, a two-hybrid 
assay was performed. In the wild-type strain, the constructs 
LexA-yGPNl, VP16-yGPN2 and pLexA-LacZ or LexA-yGPNl, 
VP16-yGPN3 and pLexA-LacZy^&K introduced. Wild-type GPN 
genes were replaced by mutated ones (yGPNl^^j^K' y^^^-^EinK 
and yGPNSgijg^). The relative (3-galactosidase activities of the 
two-hybrid assays are shown in Figure 2A. These two-hybrid 
assays clearly demonstrated that the replacement of a single resi- 
due in yGPNs, namely yGPNlg,,^' yGPN2gjj2 or yGPN3E„„, 
significantly affected yGPNl|yGPN2 and yGPNl|yGPN3 inter- 
action. While only partial loss of interaction between GPN^-^, 
and mutated GPNs was observed, almost total abolishment of 
the interaction was observed when the GIu residues were mutated 
in both GPNs (Fig. 2A). Because any change in protein amounts 
between GPN partners could alter the interpretation of two- 
hybrid results, we performed a western blot analysis showing 
that the level of mutated form yGPNl^jpj^, yGPN2gjj2K ^^id 
yGPN3£jjQ^ did not decrease compared with the wild-type form 
(Fig. 2B). In summary, we demonstrated that the E110K/E112K 
substitutions markedly alter the in vivo interaction between 
yGPNl|yGPN2 and yGPNl|yGPN3. 



The interaction between GPNs is essential for their cel- 
lular function. To establish whether the yGPNl|yGPN2 and 
yGPNl|yGPN3 interactions are important for the survival of 
yeast cells, a genetic test using the plasmid-shuffling method 
with 5-FOA selection was performed. Haploid strains, which 
lack the genomic copy ofyGPNJ, yGPN2 oi yGPN3, were rescued 
by the URA3 shuffle plasmid encoding the wild-type form of 
yGPNl, yGPN2 or yGPN3 genes, respectively. These strains were 
transformed by a second plasmid whose expression was placed 
under the control of the pGALlO galactose-inducible promoter 
and which contains appropriate WT GPNs or mutated forms 
for yGPNlgjj2> yGPN2g||2 or yGPNS^ug. Exponentially growing 
cells were spotted onto selective medium (containing 5-FOA), 
allowing to shuffle out the URA3 plasmids encoding the wild- 
type forms oi yGPNs (Fig. 3). Plated cells were grown with dif- 
ferent carbon sources such as raffinose, raffinose-galactose (98/2 
and 50/50) and galactose to trigger an increasing induction of 
GPN proteins expression. Figure 3 shows that after shuffling out 
the URA3 plasmid expressing the wild-type form of GPN, the 
copy of WT yGPNl, -2 and -3 under the control of pGALlO 
were able to complement the deleted endogenous copy for cells 
growth. However, this complementation is dependent on the 
carbon source used, as, for instance, overexpression (galactose) 
of yGPNl (Fig. 3) and underexpression (raffinose) of yGPN3 
(Table SI) drastically inhibited cell growth. Cells expressing 
yGPNlgjj2K mutant are unable to grow whatever the carbon 
source, and cells expressing yGPN2gjj2K ^^d yGPN3£j,„^ mutants 
grow only if GPN expression is strongly induced with galactose 
(50% or 100%) (Fig. 3). We noted that the cell growth rates of 
these yGPN2gjj2K ^rid yGPN3£,,(,^ mutants are nonetheless sig- 
nificantly reduced compared with the wild-type form in liquid 
culture. Another substitution (yGPN2^jj,^ and yGPN3gj|(,^) was 
also tested, as an alanine residue was predicted to have a milder 
effect on the yGPNl|yGPN2 and yGPNl|yGPN3 interactions 
than a lysine. However, similar results as those described above 
for GIu substitutions to Lys were obtained (Table SI). To address 
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Table 2. Stereochemical and energetic analysis of yGPN comparative models 



Score 


PAB0955 


yGPNI 


yGPN2 


yGPNS ^ 


PROSA 


-9.32 


-6.37 


-6.96 


-7.44 


HBOND 


398 


383 


366 


353 i 
1 


PROCHECK 


91.7 


90.5 


89.2 


93.5 


DDFIRE 


-669.1 


-604.9 


-613.7 


-598.9 ] 


FRST 


-33162.4 


-11619.5 


-14911.3 


-18082.2 


XPLOR 


-447.3 


-459.6 


-351.7 


-570.9 ] 


IVIOLPROBITY 


2.42 (52"" percentile) 


2.28 (60* percentile) 


2.44 (50* percentile) 


2.21 (64* percentile) 


Ca-RMSD 


0.0 


0.141 (241) 


0.173 (236) 


0.186(236) 1 



Values for PROSA are Z-scores obtained from the Prosa2000 package'* (lower values indicate better structures). Values for HBOND show the total num- 
ber of detected hydrogen bonds using heavy atoms of the structures (in-house tool, higher values indicate better structures). Values for PROCHECK 
indicate the percentage of dihedral angle (ip,i|)) values located within the most favored regions" (higher values indicate better structures). Values for 
DDFIRE are obtained from the DFIRE-pseudo potential function'" in units of kcal/mol (lower values indicate better structures). Values for FRST are 
obtained from the pseudo-potential energy FRST*" in pseudo units (lower values indicate better structures). Values for XPLOR are those from the non- 
bonded van der Waals potential energy terms using Charmm22 force field" in unit of kcal/mol (lower values indicate better structures). Values for Mol- 
Probity were obtained from the web server. ''The lower the score, the better the structure quality; 100* percentile represents an ideal protein structure 
whereas 0* percentile represents the worst case. Values for Ca-RMSD obtained from the sup3d global superposition™ indicate the RMSD between Ca 
atoms of the PAB0955 template and Ca atoms from each model. The number of parenthesis indicated the number of superimposed residues that have 
their Ca-RMSDs below the given value. The total alignment length is 274 residues. 



Table 3. Mutational energetic cost for Glu-^Lys replacement in 
yGPN|yGPN interactions 



AE (kcal/mol) 


yGPN1„™lyGPN2 


-1-32 


yGPNI |yGPN2,„,„ 


-1-139 


yGPN1,„„K,lyGPN2„„,„ 


4-220 


yGPN1,.„JyGPN3 


-1-98 


yGPNI |yGPN3,,„,„ 


-1-177 


yGPN1,„„JyGPN3,^„„,, 


4-225 



the specificity of the Glu mutation on cell viability, additional 
mutations were tested on residues proximal to the Glu residue. 
The highly conserved residues Gin, Ile/Val were chosen, and 
the following mutants generated: yGPNl^jj^^, yGPNl^jj^, 
yGPN2^„„^, yGPN2^,„^, yGPNS^.^,^ and yGPN3„„,^. From the 
molecular models, these mutations were predicted to have minor 
effect on the yGPNl|yGPN2 and yGPNl|yGPN3 interactions. 
Experimental results show that, indeed, these mutations do not 
affect cell viability and cell growth whatever the carbon source 
used (Table SI). In conclusion, glutamate residues at position 112 
in yGPNI and yGPN2 and 110 in yGPN3 are crucial for both 
the interactions between yGPNl|yGPN2 and yGPNl|yGPN3 
and cell viability. Otherwise, mutations in the same region 
not predicted to interfere with these interactions do not affect 
cell viability. This strongly suggests that yGPNl|yGPN2 and 
yGPNl|yGPN3 interactions are crucial for cell viability. 

yGPN3 is involved in chromosome segregation. To assess 
whether yGPN3 plays a role in sister chromatid cohesion, a 
cohesion assay with the GFP-tagged chromosome strategy as 
described earlier'* was performed on yGPN3. Because yGPN3 was 
previously found to be an essential gene,'^ we first constructed 
a conditional AyGPN3 strain by introducing in S. cerevisiae, a 
plasmid encoding the wild-type yGPN3 gene whose expression 
was placed under the control of the galactose-inducible promoter 



pGALlO and then inactivating the chromosomal copy. To GFP- 
tag the chromosome of the wild-type BY4742 and AyGPN3 iso- 
genic strains, tet0112 and TetR-GFP constructions were inserted 
at the URA3 and LEU2 loci, respectively, in both strains. 
Chromatid cohesion was examined after nocodazole treatment 
to block cells at the G,/M transition. Wild-type and yGPN3- 
regulatable mutant cells were grown to mid-log phase, exposed to 
nocodazole, and then yGPN3 higher expression was induced with 
the addition of galactose. After 30 and 60 min of induction, cells 
were harvested, and sister chromatid cohesion was estimated by 
counting among the doublet cells with a large bud, those having 
one or two GFP dots. Figure 4 shows representative microscopy 
observations where the GFP dots are clearly distinguished. The 
ratio of cells with premature separation of sister chromatids is 
plotted for both strains as a function of time (Fig. 4). For the 
wild-type strain, the ratio of cells with two distinct GFP dots was 
below 8% and remained stable after galactose induction. On the 
contrary, the mutant strain exhibited an increase in the number 
of cells with two-distinct GFP dots from 11% at t = 0 to 28% 
at t = 60 min after induction. We concluded from these data 
that higher expression of yGPN3 leads to a severe defect in sister 
chromatid cohesion, and, thus, the function of yGPN3 is related 
to this important mechanism. 

Discussion 

The GPN-loop GTPases genes are present in a single copy in 
almost all archaea and systematically in three copies (GPNl, 
GPN2 and GPN3) in all sequenced eukaryotic genomes. The 
deletion of GPNl, GPN2 or GPN3 gene in S. cerevisiae results 
in lethality.'* This points out that these three GTPases, unable to 
complement each other, perform separate functions in yeast cells. 
The analysis of a multiple sequence alignment and phylogram 
tree, constructed with GPNl, GPN2 and GPN3 genes from taxo- 
nomically distant species (Figs. S2 and 3), suggests that these 
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Figure 2. Effect of Glu-> Lys substitutions on GPNs interactions. (A) Two-hybrid assays were performed to measure the level of interaction between 
yGPNI |yGPN2 and yGPNI |yGPN3. Constructs tested are indicated in the bottom panel, and interaction data from yeast two-hybrid assays are indicated 
in relative (J-galactosidase units for each combination of GPNs. WT and mutant. E-'-K symbolize the Glu residue substitution to Lys (Ell 2 for yGPNI 
and yGPN2 and E110 for yGPN3), respectively. Error bars represent SD from the mean value of the series of five individual measurements, which were 
representative of at least two independent experiments performed on two distinct clones. (B) Western blotting (WB) of WT and mutant GPNs from 
two-hybrid assays, using anti-LexA antibody (Santa Cruz) for GPN1 detection and anti-VP16 antibody (Santa Cruz) for GPN2 and GPNS). Loading con- 
trols (LC) stained with coomassie blue SafeStain (Invitrogen) are shown (bottom panel). It shows that the expression level does not decrease in mutant 
proteins. 



GPN-loop GTPases have a common ancestor. The existence of 
functional divergence and a common ancestor^^ between GPNl, 
GPN2 and GPN3 point out that these GPN-loop GTPases are 
related by duplication and therefore are paralogous genes. 

Systematic yeast two-hybrid assays performed on yeast 
GPNs""'^ suggested that some GPNs interact with each other. 
Using TAP-tag purification technique, yGPN|yGPN interactions 
were detected in yeast: GPN1|GPN2'^ and GPN1|GPN3 (this 
work). GPN1|GPN3 complex was also identified in human. 
The need to use a cross-linker in this work to freeze the interac- 
tions indicates that the association between the GPNs is probably 
transient and/or weak. The nature of this interaction suggests 
the probable existence of interchangeable interactions of yGPNI 
with yGPN2 or yGPN3. In addition, neither two-hybrid sys- 
tems could detect GPN1|GPN1, GPN2|GPN2, GPN3|GPN3 
and GPN2|GPN3,"'"' nor could TAP-tag purification detect 
GPN2|GPN3 complexes. '^'^^ Knowing the dimeric structure of 
the GPN-loop GTPase PAB0955 from P. ahyssi," it was hypoth- 
esized that yGPNI may interact with yGPN2 or yGPN3 using 
a similar dimeric assembly. The in silico modeling revealed that 
the dimerization interface of PAB0955 was conserved in each 
yeast GPN-loop GTPase. Scrutiny of the GPNjGPN complexes 
revealed that a strictly conserved residue: yGPNlg^^) yGPN2gjj2 



or yGPN3£jj„, as well as PAB0955£,„7 ''^ template, is critical in 
the energetics of the interface interaction (Fig. 1). The structure 
analysis of GPNs reveals that this Glu residue is solvent-exposed 
in the monomeric form, and its replacement by the Lys or Ala resi- 
dues is not predicted disruptive for the subunit itself This predic- 
tion is consistent with the similar protein expression level observed 
between the WT and the Glu residue mutant forms of GPN using 
western blot assay (Fig. 2B). The yeast two-hybrid assay on vari- 
ous constructions: yGPNl^.j-|yGPN2,jj,^, yGPNl,jj,^|yGPN3^p 
yGPNl,„,JyGPN2^^, yGPNl^JyGPN2,„,, yGPNl,„,Jy 
GPN2,,,,^,yGPNl,„,JyGPN3^^, yGPNl^^|yGPN3,„„^ and 
yGPNlgjj2KlyGPN3£|jg^ revealed that the interaction is weaker 
when the Glu residue is absent from one of the two interacting 
partners and abrogated when the Glu residue is absent from both 
partners. The viability of mutants carrying the Glu substitutions 
was examined using 5-FOA shuffling assay. Glu substitution to 
Ala or Lys residue confers a lethal phenotype for yGPN, which 
can partially be rescued by overexpression in case of yGPN2 and 
yGPN3 (Fig. 3; Table SI). Finally, to investigate if the local back- 
bone conformation carrying the Glu substitution is not critical 
for the structure of yGPN proteins, conserved proximal residues 
were also mutated to Ala: yGPNl^jj^,^, yGPNl^^^, yGPN2(j||(,^, 
yGPN2yjjj^, yGPN3Qj„g^ and yGPN3,j„g^. In all these mutants. 
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Figure 3. Plasmid shuffle complementation tests of GPNs wild-type (WT) and mutants (El 12K for yGPNI and yGPN2 and El 10K for yGPN3). Cells were 
grown on control (C) or 5-fluoroorotic acid (FOA) medium. Cells were incubated for 3 d at 30°C with different carbon sources, such as raffinose, raf- 
finose/galactose (98/2 and 50/50) and galactose, to trigger an increasing induction of GPN proteins expression. 



a normal growth phenotype was observed. Taken altogether, 
these observations show for the first time that this Glu residue 
is important for the GPNs interactions and strongly suggest that 
these interactions are crucial for GPN function. It begs the ques- 
tion whether this dimeric structural assembly is unique or not 
to GPN-GTPases. The GPN GTPases belongs to the P-loop 
containing nucleoside triphosphate hydrolases superfamily (class 
of a/ (3 proteins, SCOP database^''). In this superfamily, the cell 
division regulator MinD appears in many ways similar to the 
GPNs structure. Despite a low sequence Identity between MinD 
and PAB0955 that was estimated to be around 6% according to a 
structure-based alignment made with sup3d,^° it was observed that 
the core of MinD structure^' superimposes onto that of PAB0955 
structure with a root-mean-square deviation of 1 A for 79 core Ca 
atoms (Fig. S4). MinD undergoes ATP-dependent dimerization 
in solution,^'' and the crystal structure reveals that a conserved 
Glul26 residue is located at the dimer interface. This residue is 
located right next to the G3 box (Fig. SI) and corresponds to 
Glul07 in PAB0955 (Fig. S4). In conclusion, convergent evolu- 
tionary and molecular data strongly suggest that the interactions 
of GPN-GTPases paralogs are of heterodimeric nature in the cell. 

It has been recently shown that yGPNI and yGPN2 are 
involved in sister chromatid cohesion mechanisms. A system- 
atic analysis to identify new proteins involved in maintenance 
of chromosome stability revealed that the temperature-sensitive 
yGPN2 mutant exhibited sister chromatid cohesion defects and 
was required for establishment of sister chromatid cohesion.' 
Regarding yGPNI, it was established that cells overexpressing 
this protein during anaphase displayed defects in sister chromatid 
cohesion.** Furthermore, repression of yGPNI expression strongly 
disturbs cell cycle progression with an important S-phase delay 
and an abnormal timing in nuclei migrations during mitosis.** Up 
to date, the function of yGPN3 had remained unknown. Here, 
we tested for the first time the implication of yGPN3 in this 



chromatid cohesion mechanism and discovered that, like yGPNI, 
overexpression of yGPN3 induces sister-chromatid defects. The 
fact that yeast GPN proteins have non-redundant functions that 
all participate in sister-chromatid cohesion mechanism demon- 
strates that they play a significant role in this cellular process. 

Before being described as involved in sister chromatid cohe- 
sion mechanism, GPNs were first described in stable association 
with RNA polymerase subunits (refs. 15, 17, 23, 27 and 28) and 
transcription regulators.^'-" Later, studies refined this activity and 
demonstrated that GPNl and GPN3 were essential to RNAPII 
nuclear import.'^'^^'^" Actually, with the formations gathered on 
the GPN-GTPases, it appears currently difficult to establish a 
correlation between the association of GPN-GTPases with RNA 
polymerase II and GPN-GTPases function in sister chromatid 
cohesion mechanism. However, recent works describe the exis- 
tence of intertwining between these two mechanisms. For exam- 
ple, a DNA loop, allowing physical proximity between enhancer 
and promoter for activation of gene transcription, can be stabilized 
by the cohesion complex.^' We may hypothesize that the GPN- 
loop-GTPases, in addition to their chromosome pairing activity, 
could play a role in the connection of these two mechanisms. 

Materials and Methods 

Conditions for yeast growth. For chromatid cohesion assays, 
cells were grown in yeast extract/peptone/raffinose medium, sup- 
plemented with 4% of yeast extract/peptone /galactose medium 
at 30°C. For G^/M synchronization, cells were treated with 
20 |xg nocodazole (Sigma-Aldrich) per mL of medium for 3 h 
(doubling time). For galactose induction, 5% of galactose at 
40 g/L was added. For TAP-tag experiment (Open Biosystems), 
yGPNl-interacting partners were isolated from yeast grown in 
2 L of yeast extract/peptone/dextrose liquid medium until an 
OD at 600 nm of 0.8 was reached. For two-hybrid assays, cells 
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were cultured in raffinose synthetic complete medium minus 
uracile and histidine (MP Biomedicals) supplemented with 0.2% 
of galactose synthetic complete medium minus uracile and histi- 
dine at 30°C. Yeast proteins for western blot were extracted from 
cells grown in raffinose synthetic complete medium minus his- 
tidine supplemented with 10% of galactose synthetic complete 
medium minus histidine at 30°C. For survival test on 5-FOA, 
cells were cultured in raffinose or raffinose/galactose (98/2) (for 
yGPN3 construct) synthetic complete liquid medium minus 
uracile and histidine. Cells with suitable dilution were spread on 
5-FOA 2g/L agar plates containing raffinose, raffinose/galactose 
or galactose synthetic complete medium minus histidine. Plates 
were incubated for 2-3 d at 30°C until colonies were grown. 

Construction of plasmids and yeast strains. All the strains 
and plasmids used in this study are listed in Table S2. All the 
constructed and transformed strains are derived from the 
sequenced strain S288C. For the construction of conditional 
alleles of yGPN3, a PCR-based gene deletion strategy was used to 
generate a start-to-stop codon deletion oi ylr243W. The disrupted 
gene was replaced by a hygromycin-resistance marker gene by in 
vivo homologous recombination. The deletion was performed 
into a BY4742 strain (AD5) complemented with a rescue plasmid 
(pSBTN-AE68) expressing the ORF ylr243w from the ADHl 
promoter. The hygromycin gene was amplified from the plasmid 
pAG32 using the primers oAM29 (5'-CTT AAC AAG TAC 
AAA TAG ACT AAT CAG CAT TGG AAA ATC AAT ACA 
GCT GAA GCT TGG TAC GG-3') and oAM30 (5'-CTT TTT 
ATA TGA ATC AAG CGT ACA TAA TTT TCT CTA TAA 
GCA TAG GCC ACT ACT GGA TCT G-3'). The resuking 
DNA fragment contained 40 bp terminal ends homologous to 
flanking regions containing the respective translation start and 
stop sites oi ylr243w. Cells were selected after transformation on 
plates containing rich medium supplemented with 200 |xg/mL 
of hygromycin. Three colonies from independent experiments 
were selected and shown by PGR amplification to contain the 
A)//r2-^3w::FIygroR allele correctly substituted at the ylr243w 
locus. tetR-GFP and tet0112 tagging chromosome for chroma- 
tid cohesion assays was previously described (Alonso et ak, 2011). 

Mutations into yCPNl, yGPN2 and yGPN3 (yGPNl^^^K' 
yGPNl,„,„ yGPNl„„,, yCPNl^,,,, yGPN2,,,,, yGPN2,„,„ 
yGPN2„„„ yGPN2Q,,„„ yGPN3,„„„ yGPN3,„„,, yGPN3„„,, 
and yGPN3Qj„3^) were introduced using the three fragment-homol- 
ogous recombination system.^^^ The lexA DNA-binding domain 
and VP16-activating domain of the hybrid transcription factor 
LexA-VPl6 were amplified by PGR and fused to yCPNl, yGPN2, 
yGPN3 and yCPNl^^^,^, yGPN2^„^^, yCPMS^^^^ (LexA-yGPNl, 
LexA-yGPNlg„2K VPl6-yGPN2, Y?\G-yG?m^^^^^, VP16- 
yGPN3 and VPl6-yGPN3Ej,„J. The different yGPN gene com- 
binations tested for two-hybrid interactions were co-expressed in 
the same plasmid under the control of the inducible bidirectional 
promoter pGALl-10. All constructs and mutations were checked 
by DNA sequencing. Each vector was co-transformed in yeast with 
p8op-LacZ multicopy vector (Glontech) {Ura3, 2\x). 

BY4741 isogenic derivative strain expressing yGPNl 
C-terminally tagged with a tandem affinity purification 
(TAP) -tag was obtained from Openbiosystems. 




0 

0 30 60 

time (min) after galactose addition 



Chromatid cohesion assays. Assays were performed as previ- 
ously described.^ The control strain was AI06 and the yGPN3 
mutated strain was pSBTN-AJ87. 

Comparative modeling of yeast GPN proteins. Comparative 
models (targets) were built upon the protein template of dimeric 
PAB0955 crystal structure [PDB lYRA" using a semi-automated 
protocol (unpublished)]. Targets yGPNl, yGPN2 and yGPN3 
and template were aligned using ClustalW.^^ The alignment was 
manually refined using the INTERALIGN program^'' in which 
insertions/deletions were adjusted to avoid secondary structure 
elements. In addition, the sequence alignment was adjusted so 
that all the conserved G-boxes found in GTPase families were 
aligned. Sequence identities between yGPNl, yGPN2 and 
yGPN3 with PAB0955 template are 21%, 21% and 18%, respec- 
tively. Deletion in solvent-exposed protein loops were constructed 
using sub-optimal loop alignments followed by a loop closure 
protocol.^" Sub-optimal alignment means that the exact position 
of a loop deletion is varied within a user-given range: usually 
+/- 3 residues. Short insertions (< 17 residues) were modeled using 
RAMP,'' whereas the very long insertion in yGPNl was interac- 
tively modeled as described elsewhere.^'' Cartesian positions of 
all atoms were energy-minimized as previously described."" The 
quality of comparative models was evaluated using Procheck,'' 
Prosa,''* DFIRE,'' FRST*° and MOLPROBITY.^' 

yGPNl|yGPN2 and yGPNl|yGPN3 heterodimers were mod- 
eled by superimposing yGPNl coordinates onto those of chain 
A of PAB0955, whereas yGPN2 and yGPN3 coordinates were 
superimposed onto those of chain B of PAB0955. The superposi- 
tion was performed using three conserved |3-strands of PAB0955: 
18-21, 46-50 and 112-117. Side-chains of heterodimers were 



Figure 4. Cohesion defects Increase upon yGPN3 overexpresslon. The 
graph represents the percentage of cells with double GFP dots In strain 
arrested In metaphase with nocodazole. yGPN3 expression was trig- 
gered by addition of galactose In the medium and cells were counted 
30 and 60 mIn after galactose addition. Error bars represent SD from 
the mean value of a series or three Independent experiments. In each 
experiment, counting was performed on 100 cells. 
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optimized using in-house programs,^" checked with the Xfit 
viewer/^ and the positions of atomic coordinates were refined 
using X-PLOR*^ using the all atom force field Charmm22.'''' The 
GDP molecules bound into dimeric PAB0955 were inserted into 
all yGPN models without further refinement. Buried surface area 
was computed using the buried program.'" 

The mutational energy cost for replacing Glull2 residue was 
performed by computing AAE^^^ = AEj^ - AE^,'"' where AEj^ 
and AEg are the relative potential energies of the mutant Lys 
and WT Glu in an optimized protein environment of the crystal 
structure, respectively.^' The potential energy is the sum of vdw 
and Coulombic interactions. All the relevant energy parameters 
were taken from the Charmm22 all-atom force-field library. The 
double mutation (Glu^^Lys) in yGPNl-yGPN3 heterodimer 
was computed in a similar way, except that both mutations were 
simultaneously optimized. 

P-Galactosidase activity tests. The final yeast transformants 
were analyzed for the induced (3-galactosidase activity. Fusion 
protein expression of LexA-yGPNl, VPl6-yGPN2 and VP16- 
yGPN3 was driven under the pGAL inducible promoter. For 
activity quantification, the transformants were inoculated into 
3 mL liquid medium, the cells harvested at an OD^^^^ ranging 
from 0.8-1, washed in water, centrifuged and lysed by two freeze/ 
thaw cycles. Cells were resuspended in 300 |xL of Sarcozyl buf- 
fer [60 mM Na,HPO,, 40 mM NaH,PO,, 10 mM KCL, 1 mM 
MgSO,, 0.06% (vol/vol) Sarcosyl (Fluka) and 2.6 mM DTE 
(1.4-Dithioerythritol Sigma) PH 7.0] and gently agitated 30 min 
at 30°C. Ffity (xl of samples were taken and 150 |xl of Sarcosyl 
buffer supplemented with 0.25 mM MUG (Methylumbelliferone 
(3-D-pyranoside Sigma) was added to quantify (5-galactosidase 
activity. After a 30 sec incubation time at 37°C and determination 
of the cell density at 600 nm, 4-MU fluorescence was read using 
a 360 nm excitation filter and a 460 emission filter. Fluorescence 
and absorbance were measured using a Fluostar optima (BMG) 
spectrofluorimeter. One FU/sec/OD indicated in the graph corre- 
sponds to 2.3 X 10' molecules of MU/second/cell. (3-Galactosidase 
values reported are representative of two independent transfor- 
mants tested for each construction introduced in yeast. 

Western blot. Cell disruption and protein extraction were 
performed as previously described (Alonso et al., 2001). A total of 
10 or 20 |JLg of total proteins were loaded per wells onto a 4-12% 
Bis-Tris NuPAGE precast electrophoresis gel (Invitrogen). After 
migration, the proteins were transferred onto a PVDF membrane 
(Millipore). The blot was then probed with an anti-LexA (2—12, 
Santa Cruz Biotechnology, 1/200) or anti-VP16 antibody (14-5, 
Santa Cruz Biotechnology, 1/200). Signal was developed using 
the immobilon chemiluminescence kit from Millipore. 

TAP-tagged GPNl purification, SOS-PAGE, NanoLC-MS/ 
MS analysis, MS/MS assignment and statistical procedures. 
For TAP-tagged protein purification and SDS-PAGE analysis, 
yeasts were grown in 2 L of YPD liquid medium until an OD 
at 600 nm of 0.8 was reached. For in vivo formaldehyde cross- 
linking, two different concentrations of formaldehyde (0.15% 
and 0.2%) were directly added to the yeast cell culture, and 
cells were incubated at 30°C for 10 min. The cross-linking reac- 
tion was quenched for 20 min at 30°C by glycine with a final 



concentration of 0.125 M. After cross-linking, cells were col- 
lected and washed three times with 100 mL of ice-cold water. 
Pellets were frozen in liquid nitrogen and stored at -80°C until 
use. Cell pellets were lysed mechanically by means of steel beads 
using a mixer mill from Retsch in cold conditions. The purifica- 
tion of the TAP-tagged protein and its partners was done starting 
from cell lysates containing about 40 |JLg/ml total proteins. The 
purification was performed as described by Rigaut et al.**^ in the 
presence of benzonase (25 U/mL) and without EDTA. For the 
final elution step, proteins retained on the calmodulin beads were 
eluted by addition of 70 |jlL of lithium-dodecyl-sulfate solution 
(Invitrogen). After 5 min incubation at 95°C, the whole sample 
was resolved by SDS-PAGE onto a 4-12% bis-TRIS NuPage gel 
(Invitrogen) with 3-(N-morpholino) propanesulfonic acid buffer 
(Invitrogen) as running buffer, but with a short migration time. 
The whole protein content from each lane was excised as a single 
polyacrylamide band, destained with water, treated and proteo- 
lyzed with trypsin as described previously.'"* 

For nanoLC-MS/MS analysis, NanoLC-MS/MS experiments 
were performed on a LTQ-Orbitrap XL hybrid mass spectrometer 
(ThermoFisher) coupled to an UltiMate 3000 nano-LC system 
(Dionex-LC Packings) essentially as described in reference 49. 

For MS /MS assignment, peak lists were generated with 
the MASCOT DAEMON software (version 2.2.2) from 
Matrix Science using the extract_msn.exe data import filter 
(ThermoFisher) from the Xcalibur FT package (version 2.0.7) 
from ThermoFisher. Data import filter options were set at: 
400 (minimum mass), 5,000 (maximum mass), 0 (grouping 
tolerance), 0 (intermediate scans) and 1,000 (threshold). 
MS/MS assignments were performed using the MASCOT search 
engine (version 2.2.2, Matrix Science) against a local database 
constructed with the S. cerevisiae protein sequences from the 
Saccharomyces Genome Database (SGD release 20100105) avail- 
able at www.downloads.yeastgenome.org. This database release 
comprises 5,885 protein sequences, totaling 2,916,234 amino 
acids. Searches for trypsic peptides were performed with the fol- 
lowing parameters: a mass tolerance of 5 ppm on the parent ion 
and 0.5 Da on the MS/MS, static modifications of carbamido- 
methylated Cys (+57.0215) and dynamic modification of oxidized 
Met (+15.9949). The maximum number of missed cleavages for 
trypsin was set at 2. All peptide matches with a peptide score 
above its peptidic identity threshold set at p < 0.05 and rank 1 
were filtered by the IRMa 1.23 software.'" 
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